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PREFACE 


The  work  presented  in  this  report  has  been  sponsored  by  the 
Federal  Aviation  Administration  as  part  of  a program  to  provide 
navigational  aids  for  the  safe  landing  of  aircraft.  The  program 
has  been  concerned  with  the  application  of  instrument  landing 
aids,  and  in  particular  the  development  of  models  to  predict  the 
performance  of  localizer  and  glide  slope  antenna  systems.  The 
present  report  is  concerned  with  the  glide  slope  portion  of 
instrument  landing  system,  and  the  effects  of  airport  topography 
on  its  performance. 
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1.  INTRODUCTION 


A mathematical  model  has  been  developed  for  predicting  the 
performance  of  image-type  glide  slope  arrays.  The  basic  theory  is 
developed  in  the  following  section  of  this  report.  In  Part  A of 
that  section  a model  of  the  glide  slope  multipath  problem  is 
developed  while  in  Part  B the  techniques  for  dealing  with  irregular 
terrain  profiles  are  outlined.  In  the  final  section  of  the  report 
some  illustrative  numerical  results  are  presented. 


2.  BASIC  THEORY 
PART  A 

The  presence  in  the  airport  environment  of  such  large  man- 
made structures  as  aircraft  hangars  as  well  as  such  natural  terrain 
features  as  hillsides  can  lead  to  glide  slope  course  derogation 
through  multipath  scattering.  We  begin  our  treatment  of  glide 
slope  siting  problems  by  developing  a model  for  predicting  the 
amount  of  such  multipath  derogation.  A typical  situation  is 
depicted  in  Figure  1. 

The  t uchdown  point  on  the  centerline  of  the  runway  opposite 
the  glide  slope  array  is  chosen  as  the  origin  of  coordinates  0. 

The  z-a>is  is  chosen  to  be  the  vertical  axis  passing  through  the 
origin  0,  while  vhe  x-and  y-sxes  are  parallel  to  and  perpendicular 
to  the  centerline,  respectively.  The  positive  z-axis  points  out  of 
the  page  (Figure  1).  The  ground  plane  is  assumed  to  be  perfectly 
conducting.  Consequently,  any  deviations  from  nominal  glide  slope 
performance  are  attributable  to  the  multipath  scattering  produced 
by  the  various  structures  (both  natural  and  man-made).  For  simpli- 
city, we  will  assume  that  all  such  structures  are  perfectly  con- 
ducting. 

In  the  situation  depicted  in  Figure  1,  the  scatterer  (perhaps 
a small  hill)  is  illuminated  by  the  glide  slope  array  and  scatters 
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FIGURE  1.  TYPICAL  GEOMETRY  FOR  MULTIPATH  SCATTERING 
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some  of  this  signal  back  into  the  path  of  the  approaching  aircraft. 
We  now  proceed  to  calculate  mathematical  expressions  for  these 
scattered  fields  at  the  receiver.  For  convenience,  we  will  work 
with  magnetic  field  intensities. 

For  a perfect  conductor,  the  surface  current  density  K is 
given  by: 

T * n x if  , (1) 

where  n is  the  local  unit  normal  vector  pointing  out  of  the 
scattering  surface  and  H is  the  total  magnetic  field, 

TT  - IT.  ♦ Hs,  (2) 

where  H.  is  the  incident  and  H the  scattered  field  on  the  surface 
1 s 

of  the  scattering  structure.  In  terms  of  the  surface  current  den- 
sity K,  the  scattered  field  Hg  at  the  receiver  is  given  by  the 
following  surface  integral: 

Kfy  <+  /CTfT^Tj.TOd,  (3) 

s 

In  Equation  (3),  the  vectors  r^  and  r denote,  respectively,  the 
position  vector  of  the  receiver  and  the  position  vect  ,r  of  some 
arbitrary  source  point  cn  the  surface  of  the  scatterer  relative  to 
the  origin  of  coordinates  0.  The  two-point  Green's  function  G("rj,"r) 
is  given  by: 

Gffj.7)  - (4) 

where  k » 2 n/ A and  X.  is  the  wavelength  of  the  incident  radiation. 

The  integral  is  taken  over  the  surface  S of  the  scatterer. 

We  will  adopt  here  the  single  scattering,  physical  optics 
approximation  for  the  current  density  K.  Specifically,  we  will 
assume  that  on  those  portions  of  the  scattei  ng  surface  not  directly 
illuminated  by  the  glide  slope  array,  K is  identically  zero  and 
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that  on  the  directly  illuminated  portions  of  the  scattering  surface, 

K is  proportional  to  twice  the  tangential  component  of  the  incident 
magnetic  field: 

F = 0 on  S_  (5) 

F - 2 (n  x TT)  on  S+  '6) 

where  and  S_  denote,  respectively,  the  illuminated  and  unillu- 
minated portions  of  the  scattering  surface  S (see  Figure  1).  For 
the  current  distribution  defined  by  Equations  (5)  and  (6),  Equation 
(3)  becomes. 

lTs(ri)  n (7 ) xH*  (7) ^ x 7 G(‘rJl7r)ds.  (7) 

S+ 

The  physical  optics  scattering  model  represented  by  equation  (7) 
is  based  upon  the  existence  of  a sharp  shadow  boundary  on  the 
surface  of  the  scattering  structure.  This  assumption  ignores  such 
phenomena  multiple  reflections,  surface  waves,  and  diffraction  but 
should  provide  reasonably  accurate  results  as  long  as  the  surface 
features  of  the  scattering  structure  do  not  vary  greatly  over 
dimensions  which  are  small  compared  to  the  wavelength  A (about 
3 feet  at  glide  slope  frequencies). 

In  principle,  the  solution  to  the  problem  at  hand  requires 
only  that  the  magnet  field  be  determined.  For  simplicity,  we 
will  assume  that  the  glide  slope  array  is  made  up  of  electrically 
short  dipoles.  With  this  assumption,  we  in  effect  approximate  by 
a cosine  distribution  the  actual  azimuthal  pattern  of  the  half-wave 
dipoles  which  make  up  glide  slope  arrays.  This  approximation  should 
not  be  too  restrictive  and  it  docs  considerably  expedite  our  calcula- 
tions. Let  (0,  y , h)  denote  the  x,  y,  and  z coordinates  respectively 

A 

of  a typical  dipole  in  the  array  and  let  (x,  y,  z)  denote  the  coor- 
dinates of  a point  r at  which  we  wish  to  know  the  magnetic  field 
intensity  Ha  produced  by  the  dipole  at  (0,  ya,  h) . In  our  earlier 
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glide  slope  performance  prediction  report,*  it  is  shown  that  H (r) 

3 

is  given  by 


(8) 


where  e , e , and  e denote  unit  vectors  in  the  x,  y,  and  z direc- 
X y z 

tions  respectively,  JQ  is  a parameter  which  measures  relative  phase 
and  amplitude,  and  Dj  is  the  distance  from  the  antenna  to  the  obser- 
vation point  r: 


Dj  -C*2  ♦ (yya)2  * (z^O1* 


(9) 


To  take  into  account  reflections  from  the  ground  plane  (z  * 0) 
we  use  simple  image  theory  to  obtain  the  image  field  H(r): 

3 


«aCr) 


ikJ0 

4ir 


kDj 

TT-  O*  * (2+h)ex  J 
U2 


(10) 


where  D2  is  the  distance  from  the  image  of  the  transmitting  dipole 
to  the  field  point  r: 


°2  -Cx2  + (y-ya)2  * (z+h)234 


(11) 


The  total  field  intensity  H.  at  r is  just  the  sum  of  the  direct 
(H„)  and  ground  reflected  (H  ) signals: 

3 3 

Hi (7)  - V)  +7*(7)  (12) 

Note  that  on  the  ground  plane  (z  = 0),  the  z component  of  H.  vanishes 
identically. 
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The  field  H.  determined  by  Equations  (8),  (10),  and  (12)  can 
now  be  substituted  into  Equation  CM  and  the  integral  can  be  per- 
formed over  the  illuminated  surface  of  the  scattering  structure  to 
give  the  scattered  field  at  the  receiver  location  Tj.  Actually,  only 
the  z component  of  the  scattered  field  need  be  calculated  since  the 
receiving  antenna  responds  primarily  to  the  .sorizontal  component  of 
electric  field  which  is  proportional  to  the  z component  of  the  magnetic 
field.  It  should  be  noted  that  for  a complete  solution,  the  contri- 
bution of  the  image  of  the  scattering  structure  must  also  be  cal- 
culated. We  will  assume  that  the  receiver  is  always  in  the  far 
field  of  the  scattering  structure  so  that  the  following  asymptotic 
form  of  the  gradient  of  the  Green's  function  is  applicable: 

7ti(rj.  , * - ik  (rj-r)  e>kirl'r! 


Substituting  Euqations  (8),  (10),  (12)  and  (13)  into  Equation 
(7)  and  taking  the  z-component  of  the  resulting  vector  equation, 
we  obtain  the  following  expressions  for  Hsj(rj),  the  z-component  of 
scattered  magnetic  field  at  the  receiver  point  r{ : 


— k J 


Hsz<V 


8 tr 


c«. 


(14) 


where 


>■/; A 


,ik(VDl)  d5 


S.D1  R1 


(15) 


f F2 

(tv 


(16) 


o 


eik(R2*D1) 


ds 


(17) 


f 1 

JT7^ 

S * U1  k2 


U-f-Xr  d5  <>’> 


D1 

* 

Cx2  - (y-Yj)2  ♦ C*-h)2J  ** 

(19) 

°2 

m 

Cx2  ♦ (Y-Ya)2  * (**h)0% 

(20) 

R1 

M 

C(Xx-x)2  ♦ OTj-y)2  ♦ (z? -z)2D ** 

(21) 

R2 

m 

C(XA'x)2  ♦ (Yj-y)2  ♦ CZ,^z)2D  ** 

(22) 

F1 

9 

RjMij-y)  < QijX  ♦ n3(Z'h)  ] (Xj-x) , 

(23) 

F2 

m 

n2x  CYi-y)  ♦ QijX*n3(z*h)])  (Xj*x) 

(24) 

In  Equations  (14)  through  (24),  the  coordinates  (Xj  Yj , Zj)  denote 
the  coordinates  oi  the  aircraft  (x,  y,  z)  the  coordinates  of  some 
source  point  on  the  surface  of  the  scattering  structure  and  (nj, 
n2'  n3^  t**e  components  of  the  unit  outward  normal  to  the  surface  of 
the  scatterer  at  the  point  (x,  y,  z).  The  integrals  I3  and  1^ 
represent  the  contributions  of  the  image  of  the  scattering  structure 
in  the  ground  plane  z • 0.  Note  the  following  relationships  among 
the  four  integrals: 


*1 

(Xj,  Yj, 

zl> 

- -Ij  (Xj,  Yj , -Zj) 

*2 

>• 

X 

zl> 

- -I4  (Xj,  Yj,  -Zj) 

(25) 

Consequently,  when 

Z1  • 

0,  Hsz(rj)  = 0 as  it  should  be. 
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The  direct  numerical  integration  of  the  integrals  Ij,  , Ij, 
and  I4  for  an  arbitrary  scattering  structure  would  in  general 
require  an  inordinate  amount  of  computer  tiae.  However,  certain 
simplifying  assuaptions  can  be  aade  which  considerably  expedite  the 
evaluation  of  these  integrals.  Specifically,  we  will  assume  that 
any  scattering  surface  can  be  adequately  represented  by  a series  of 
interconnecting  plane  facets.  By  definition,  the  unit  noraal  vector 
n fill  be  constant  over  each  facet.  The  exact  size,  shape,  and 
orientation  of  these  facets  will  depend  upon  the  surface  charac- 
teristics of  the  particular  scatterer.  In  turn,  each  facet  can  be 
broken  up  into  a series  of  interconnected  rectangular  pieces.  If 
these  rectangular  pieces  are  aade  saall  enough,  the  integrals  Ij, 

I2t  Ij  and  1^  can  be  evaluated  analytically  for  each  piece  and 
then  the  results  suaaed  over  the  whole  surface  of  the  scatterer  to 
give  H$z.  To  illustrate  these  procedures,  we  will  now  evaluate 
the  integral  Ij  for  an  arbitrarily  oriented  rectangular  plate. 

In  the  course  of  this  evaluation,  we  will  discuss  the  size  res- 
trictions on  such  plates. 

In  Figure  2,  a rectangular  facet  typical  of  those  Making  up 

the  surface  of  >o ae  scatterer  has  been  drawn.  We  will  evaluate  the 

integral  I,  of  Equation  (15)  for  such  a facet  m the  Fraunhofer 
* A A 

approximation.  The  unit  vectors  n and  £ lie  in  the  plane  of  the 
facet,  are  orthonornal,  and  define  the  direction  of  the  outward 
nornal  n: 

AAA 

n • nx  ( 

The  center  point  of  the  facet  (x0,  yQ,  zQ)  will  be  used  as  a local 
origin  of  coordinates  for  the  i »tegrations  which  are  to  be  per- 
formed. The  integral  Ij  is  given  by 


I 


1 


J d/r 


2 

1 


_ik(Di*Rj) 

* 1 1 ds 


(26) 
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where 


Fj  * n2  x (Yj-y)  ♦QnjX  ♦ n3(z*h)3  (Xj-x)  (27) 

Dj  - rx2  ♦ (y-ya)2  ♦ (z-h)2]  **  (28) 

R1  * C(X1'*)2  * CYi*y>2  " (Zi-02D%  C29) 

Clearly,  the  (Hstances  Dj  and  Rj  denote,  respectively,  the  length 
of  the  vector  Dj  fro*  *he  antenna  to  a point  (x,  y,  z)  on  the  sur- 
face of  the  facet  and  the  length  of  the  vector  Rj  fro*  the  point 
(x,  y,  z)  to  the  receiver.  Fro*  Figure  2,  it  is  seen  that  tK** 
vector  Dj  can  be  represented  as  follows: 

D1  * D10  + r (30) 


where 

D10  : Xo*x  * <y0'/a)«y4(zo'h)®z 

and 

r ■ nn  ♦ ££ 


(31) 


(32) 


and  the  ranges  of  the  variables  n and  C are  - y—  < n < y-^  and 

a r At  L m m L 


- §1<  C < ft  . 


Fro*  Equation  (30),  we  have 

.2.™  1/2 


'10 


' D10  ' 


(33) 


Normally,  Djq  is  cn  the  order  of  thousands  of  feet  and  r is  of  the 
order  of  tens  of  feet  so  that  the  radical  in  Equation  (33)  can  be 
expanded  in  powers  of  r/Dj0  to  give  Dj  approximately  as  a function 
of n and  £.  To  terms  of  second  order  in  r/Djg,  we  obtain 


Dl*  D10 


♦ n cos  a *1  cos  £ ♦ yS 

zu10 


<D10 

*®10 


•r , 


(34) 


10 


where 


610  ■ ®10/D10 

(35) 

A 

cosa  * Djg’n 

A A 

(36) 

coi 6 * Djq* £ 

(37) 

The  numbers  coso  and  cosB  are  the  direction  cosines  of 
relative  to  n and  £.  In  the  Fraunhofer  approximation  which  we 
will  use,  it  is  assumed  that  the  quatratic  terns  of  Equation  (33) 
are  snail  coapared  to  a wavelength  for  any  point  (x,  y,  z)  on  the 
surface  of  the  facet  so  that  we  can  write  to  a good  approximation: 

Dj  ■ Dj0  ♦ n cosa  ♦ C cosB  (38) 

It  can  be  shown  that  the  largest  value  attainable  by  the 
quadratic  terns  in  Equation  (33)  on  the  facet's  surface  is 

A - sin2a  ♦ sinZB  ♦ ^^|cosa  cosb|  (39) 

Equation  (38)  should  provide  a good  approxination  for  the  distance 
Dj  as  a function  of  n and  £ as  long  as  An  and  A£  are  chosen  so  that 
A in  Equation  (39)  is  snail  conpared  to  the  wavelength  X.  We  will 
return  to  this  important  requirement  later. 

Similarly,  the  vector  Aj  can  be  represented  as  follows: 

#1  * A10  - r (40) 

where  AJ0  - (X1-x0)ex*(Y1-y0)ey^(Z1-z0)ez  (41) 

Thus,  in  the  Fraunhofer  approximation,  we  find 

Rj  * Rjq  - n cos  a * 5 cos6  C42) 
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A A 


COS  3 

R10  i) 

(43) 

A A 

( H i 

COS  8 

R10  £ 

where 

*10  " 

R10/R10 

(45) 

In  Equation  (42),  we  have  ignored  terms  of  second  order  and  higher 
in  r/Rjg. 


One  final  approximation  which  we  will  make  is  thatAri  and  AC 
are  sufficiently  small  compared  with  Djq  and  P^q.  that  the  variation 
of  the  amplitude  function  Fj/D^R^  in  Equation  (26)  can  be  ignored 
This  is  normally  an  excellent  approximation.  Consequently,  we  will 
assume  that  for  any  point  on  the  surface  of  the  facet  we  can  write 

f,/",2",2  F,o/Din2Bio2  <46’ 

where  F1Q,  D1Q,  and  RjQ,  are  the  values  of  Fj , 1^ , and  Rj  at  the 
center  point  (xQ,  yQ , zq)  . 

Substituting  Equations  (38),  (42),  and  (46)  into  Equation 
(26)  yields  the  following  approximate  expression  for  1^; 

*AC/2 

Ij  = eikn(cosa  'cosy)dr)f ei^(cos£-cos8)d£ 

nioRio  -AV2  'A  C/2  f 4 7 1 


Performing  the  integrals  indicated  in  Equation  (47),  we  finally 
obtain  the  following  approximate  expression  for  I j 


4F, 


10 


eik(D10*R10> 


1 

J10  10 


(48) 


sin^  ^A*)  (ccso-cosy)^  s in  ^ k t >} 


(cos  0 -cos>  ) (cosfl-cos& ) 


1 2 


A siailar  treatment  of  the  integrals  I2,  Ij,  and  14  ultiaately  lead 
to  the  following  approxiaate  expression  for  the  scattered  magnetic 
field  intensity  at  the  receiver  due  to  the  rectangular  facet: 


sz 


J F e*k^B10*R10^ 

2*2  »ioRio  A B 

J0  F20elk(D20*R105 

sin(kAlA*)/2)sin(kBlA(/2) 

2 R2 

1 d20r10 

A1B1 

Jo  Fl0*ik(Dl0+R20) 

sin(kA2AiJ/2)sin(kB2A^2) 

***  D?o*io 

A2B2 

J F,nelk{D20*R20) 
^ o zo 

sin(kAjAf7/2)sin(kBjA^2) 

(49) 

777  JT 
2 D20R20 

A3B3 

F10*  F20 * D10'  °20’  R10  and  R20  denote  the  values 
, Rj,  and  R2  at  (xQ,  yQ,  zQ)  and  the  paraaeters  A 

B^,  Bp  and  Bj  are  defined  as  follows: 

of  Fp  F2 

* Aj ) A2 , 

A * cosa  - cosy 

B - cosB  - cos6 

(50) 

Aj  * COS&j  * co«y 

Bj  ■ cosBj  - cos6 

(51) 

A2  * cosa  - cos y j 

B2  • cosB  - cosSj 

(52) 

Aj  ■ cosaj  - cosYj 

Bj  - cosBj  - cos6j 

(53) 

and  the  angles  appearing  in  Equations  (50)  through  (53)  are  defined 
as  follows: 


AAA 

Vi  * W'y  * <Vh)ez 

<w*x  * <v>\A  * <W'z 


(54) 

(55) 
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°20  " V*  * ‘WV'V’K 


*20  ' (V\>,V(Yr>'o>V(Zl*20)<’z 


and 


, __  A 

„ D10  o °10  ^ 

COS  a * -p ttf  , cos  0 * -p 


10 


'10 


ccs  a 


-*■  A 

D:o-1? 


D 


20 

A 


cos  a * I 


l10 


, co*01 


»20* 


D 


20 


cos 


8 ■ 


— A 

R10’* 

R10 


a:id  cos  a j * 


R20 


cos 


8i  m 


— A 

R20  t 


R 


20 


(56) 

(57) 

(58) 

(59) 

(60) 
(61) 


The  total  scattered  field  at  the  receiver  is  obtained  by  summing  up 
the  contributions  from  all  of  the  rectangular  facets  which  make  up 
the  surface  of  the  scattering  structure.  Segmenting  the  surface 
of  a scatterer  into  rectangular  facets  and  then  applying  the  approxi- 
mate closed  form  solution  in  Equation  (49)  to  each  facet  represents 
a far  more  efficient  procedure  for  calculating  the  scattered  field 
at  the  receiver  than  attempting  to  directly  numerically  integrate 
Equations  (14)  through  (18). 

We  now  return  to  the  question  raised  earlier  of  the  size  res- 
trictions upon  the  facets.  Clearly,  one  would  like  to  maximize 
the  area  of  each  facet  in  order  to  minimize  computer  running  time. 

The  principal  restriction  upon  the  size  of  the  facets  is  that  \nand 
AC  must  not  be  so  large  that  the  Fraunhofer  approximation,  upon  which 
the  approximate  formula  (49)  is  based,  is  violated.  That  is,  the 
quadratic  path  length  difference  terms  like  those  in  Equation  (34) 
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■ust  be  >mall  compared  to  the  wavelength  A for  all  points  on  the 
surface  of  the  facet  to  wrich  Equation  (49)  is  to  be  applied.  Thus, 
in  choo  ing  appropriate  values  for  A£  and  A£  , we  are  essentially 
faced  kith  an  optimization  problem  with  a constraint.  Specifically, 
we  want  to  maximize  t*»e  area  AnA£  of  the  facet  while  keeping  the 
quadratic  path  length  difference  terms  less  than  some  prescribed 
value.  For  example,  for  the  path  length  0^ , the  maximum  value  of 
the  quadratic  pati’  length  difference  term  is  given  by 


4 <A»?)2  2 (AO2  A17A* 

A m sin  a+ * 


8 0 


10 


8 D 


10 


COSOCOS^j 


(62) 


Depending  upon  the  sign  of  cosacosB,  this  value  will  be  achieved  at 
two  of  the  four  corners  of  the  rectangular  facet  At  all  other 
points,  the  quadratic  path  length  terms  will  be  less  than  A,  but 
always  positive.  We  can  optimize  our  choice  ofAn  andASby  maximi- 
zingA o AS  subject  to  the  condition  that  A in  Equation  (62)  be  some 
fractional  part  of  a wavelength.  This  problem  has  been  solved  in 
closed  form  using  Lagrange  multipliers.  The  optimization  procedure 
must  be  applied  to  the  quadratic  path  length  difference  terms 
affecting  all  four  distances  (Dj,  D2,  Rj , R2)  and  then  from  the 
resulting  set  of  solutions  the  smallest  values  of  An  and  are 
chosen.  . subroutine  has  been  written  to  carry  out  the  optimization 
procedure. 

Of  course,  the  scattered  fields  must  be  calculated  for  each 
dipole  in  the  array,  the  parameter  JQ  being  varied  to  take  into 
account  the  exact  phase  and  amplitude  ratios  among  the  various 
elements  and,  along  with  JQ,  the  height  h and  y-offset  ya  of  each 
element.  Our  program  automatically  offsets  the  elements  of  an 
array  to  correct  for  the  effects  of  proximity  phase  lag.  The 
y-coordinates  cf  the  elements  are  adjusted  so  that  all  dipoles  are 
at  the  same  slant  distance  from  the  touchdown  point  on  the  runway 
directly  opposite  the  array. 


IS 


t 
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Once  all  of  the  scattered  fields  have  been  calculated,  the 
calculation  of  the  difference  in  depth  of  modulation  (D.D.M)  at  a 
given  receiver  location  is  straight  forward.  The  scattered  fields 
at  the  carrier  frequency,  90Hz  modulated  frequency,  and  ISOHz 
modulated  frequency  are  first  added  to  the  direct  fields  at  those 
frequencies.  By  direct  fields  is  meant  the  fields  which  would  exist 
if  the  scatterers  were  not  present.  For  example,  if  the  ground 
plane  is  perfectly  flat,  the  direct  field  foi  any  one  of  the  three 
frequencies  at  the  receiver  location  rj  is  just  the  z component  of 
Equation  (12). 


H 


it 


(63) 


where,  here,  Dj  and  D2  denote,  respectively,  the  distance  from  the 
antenna  to  the  aircraft  and  the  distance  from  the  image  antenna  to 
the  aircraft,  and  Xj  is  the  x coordinate  of  the  aircraft.  If  the 
ground  plane  is  not  perfectly  flat  or  at  least  nearly  so  the  direct 
fields  are  calculated  differently.  These  alternate  procedures  will 
be  discussed  in  the  next  section.  In  any  case,  once  the  total  com- 
plex field  amplitudes  (direct  plus  scattered  summed  for  all  dipoles 
in  the  array)  have  been  found,  the  DDM  is  calculated  as  follows: 
let  Hc,  H1S0,  and  H90  denote  the  total  complex  amplitudes  at  the 
carrier  frequency,  thr  150  Hz  modulated  frequency,  and  the  90  Hz 
modulated  frequency.  In  terms  of  these  complex  amplitudes,  the 
difference  in  depth  of  modulation  is  given  by: 


D.D.M. 


Re 


(64) 


where  Re  denotes  the  'eal  part  of  the  complex  number.  Full  scale 
deflection  (ISO  microamps)  corresponds  to  a D.D.M.  of  .175  for  the 
glide  slope.  Consequently,  the  C.D.I.  (course  deviation  indication) 
is  given  by 

i 

150  »u  i_. (65) 


C.D.I.  - 857.14  Re 


icroamps 
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PART  B 


In  Part  A of  this  section,  we  developed  techniques  for  esti- 
nating  glide  slope  performance  degradation  due  to  multipath 
scattering.  Once  the  scattered  fields  produced  by  the  various 
scattering  structures  in  the  airport  environment  have  been  calcu- 
lated, they  are  added  to  the  direct  fields  at  the  corresponding 
frequencies  and  from  these  total  fields,  the  C.D.I.  is  calculated 
using  Equation  f*c).  The  direct  field  at  the  particular  frequency 
is  just  the  field  which  would  be  detected  at  the  receiver  in  the 
absence  of  the  scattering  structures.  If  the  ground  plane  is  other- 
wise perfectly  flat  or  nearly  so,  the  direct  field  is  given  accord- 
ing to  simple  image  theory  by  Equation  (63). 


(63) 


where  Dj  and  D2  are,  respectively,  the  distance  between  the  receiver 
and  the  transmitting  dipole  and  between  the  receiver  and  the  image 
of  the  dipole.  On  the  other  hand,  if  the  terrain  profile  exhibits 
significant  irregularities  such  as  downgrade;  or  drop  offs  along  the 
line  of  sight  between  the  aircraft  and  the  antenna  mast,  Equation 
(63)  cannot  be  used  and  the  techniques  developed  in  our  earlier 
glide  slope  performance  prediction  model  must  be  used. 


In  the  report  ILS  Glide  Slope  Performance  Prediction,  Vol.  B 
(FAA-RD-74-157.B) , we  developed  mathematical  procedures  for  calcu- 
lating the  ground  reflected  signals  resulting  from  irregular  terrain 
profiles  such  as  the  one  depicted  in  Figure  3.  The  detailed  calcu- 
lations will  not  be  repeated  here,  However,  the  updated  computer 
program  contains  the  subroutines  for  performing  these  calculations. 
One  simply  inputs  the  shape  of  the  terrain  profile  along  the  line  of 
sight  between  the  aircraft  and  the  antenna  mast,  and  the  subroutines 
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FIGURE  3.  TYPICAL  PROFILE  OF  ONE -DIMENSIONAL  TERRAIN  VARIATION 


will  generate  the  ground  reflected  signals  at  the  carrier  frequency. 
90Hz  modulated  frequency,  and  the  150Hz  modulated  frequency.  These 
signals  are  then  combined  with  the  incident  signals  directly  from 
the  antennas  to  give  what  we  have  been  calling  the  direct  fields, 
the  fields  which  would  be  measured  at  the  reeeiver  assuming  that  no 
scatterers  (hangars,  hillsides,  etc.)  were  present.  These  direct 
signals  art  then  combined  with  the  scattered  signals  calculated 
using  the  techniques  of  Part  A and  the  resulting  total  signals  are 
used  in  Equation  (65)  to  calculate  C.D.I.  Again,  if  the  terrain 
profile  is  perfectly  flat  or  nearly  so,  Equation  (63)  can  be  used 
to  calculate  the  direct  signals. 


3.  NUMERICAL  RESULTS 

Before  presenting  some  sample  numerical  results  illustrating 
the  predictions  of  our  model,  we  will  first  review  the  characteris- 
tics of  the  three  main  types. of  image  glide  slope  arrays;  the  null 
reference,  the  sideband  reference,  and  the  capture  effect  arrays. 

The  null  reference  array  is  the  simplest  of  the  image  glide 
slope  systems.  It  consists  of  two  transmitting  antennas  whose 
heights  are  in  the  ratio  of  2:1.  The  upper  antenna  is  fed  side- 
band only  signal,  the  150  Hz  and  the  90  Hz  components  being  of 
equal  amplitude  but  180°  out  of  phase.  The  lower  antenna  is  fed 
both  carrier  and  sideband.  The  carrier  signal  is  nominally  40t 
modulated.  The  ratio  of  the  sideband  currents  in  the  upper  antenna 
to  the  sideband  currents  in  the  lower  antenna  is  typically  0.3. 

The  sideband  ratio  of  0.3  would  nominally  yield  a 1.4°  course 
width  (full  deflection  0.7°  above  or  below  the  glide  path).  Let 
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Ic.  ^i5o»  an<*  ^90  denote  the  current  amplitudes  fed  to  the  antennas 
at  the  carrier  frequency,  the  150  Hz  modulated  frequency,  and  the 
90  Hz  modulated  frequency,  respectively.  The  null  reference  array 
parameters  discussed  above  can  be  summarized  as  follows: 


Carrier  Antenna 


Sideband  Antenna 


Height  * h 


I. 


150 


90 


1 

« 0.4 
» 0.4 


Height  * 2h 


0 


150 


90 


- 0.12 
—0.12 


Note  that  all  current  amplitudes  have  been  normalized  relative  to 
the  carrier  current  amplitude  (Ic  - 1). 

The  sideband  reference  array  employs  two  transmitting  dipoles 
whose  heights  are  in  the  ratio  of  3:1.  If  the  lower  ntenna  is 
positioned  at  1/2  the  height  of  the  lower  antenna  of  the  null 
reference  array  and  if  the  upper  antenna  is  positioned  at  3/4  of 
the  height  of  the  upper  antenna  of  the  null  reference  array,  the 
same  glide  angle  is  produced.  Modulated  carrier  (401  modulated) 
is  fed  to  the  lower  antenna.  Both  antennas  are  fed  separate 
sideband  signals.  The  separate  sideband  signals  fed  to  the  two 
antennas  are  equal  in  amplitude  but  are  180°  out  of  phase.  The 
amplitude  ratio  of  the  separate  sideband  signal  to  the  carrier 
sideband  signal  is  typically  0.3.  This  ratio  produces  a nominal 
course  width  of  1.4°  as  in  the  case  of  the  null  reference  array. 
The  sideband  reference  parameters  are  summarized  below. 


Lower  Antenna 


Height 
I-  - 1 


hso  " °*28 


90 


0.52 


Upper  Antenna 
Height  =■  3h 


!c  - 0 


150 


■90 


0.12 

-0.12 
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Note  that  the  1^^  and  IgQ  current  amplitudes  given  above  for  the 
lower  antenna  represent  the  sums  of  the  carrier  sideband  and 
separate  sideband  signsls  (Ijjq  e .4  - .12  * .28,  I90  * .4  ♦ .12 
3 . S2).  All  current  amplitudes  have  been  normalized  relative  to 
the  carrier  signal  amplitude. 

The  capture  effect  glide  slope  array  consists  of  three  trans- 
mitting antennas  whose  heights  are  in  the  ratios  of  1:2:3.  If  the 
lower  and  middle  antennas  are  set  at  the  same  heights  as  the  null 
reference  antennas,  the  same  glide  angle  is  produced.  We  will  not 
treat  the  clearance  signal  which  provides  a s;** fly  up  signal  at 
low  approach  angles  but  has  little  effect  upon  the  glide  ar.gle  and 
course  width.  Concerning  the  primary  signal,  the  modulated  carrier 
is  id  to  both  the  lower  and  middle  antennas.  The  modulated  carrier 
fed  to  the  middle  antenna  has  half  the  amplitude  and  is  180°  out 
of  phase  with  the  modulated  cirrier  fed  to  the  lower  antenna.  The 
carrier  signals  are  nominally  40t  modulated.  In  addition,  all  three 
antennas  are  fed  separate  sideband  signals.  The  separate  sideband 
signals  fed  to  the  lower  and  upper  antennas  have  half  the  amplitude 
and  are  180°  out  of  phase  with  the  separate  sideband  signal  fed  to 
the  middle  antenna.  The  ratio  of  the  separate  sideband  signal  fed 
to  the  middle  antenna  to  the  carrier  sideband  signal  in  the  lower 
antenna  is  typically  0.3.  This  ratio  yields  nominally  a course 
width  of  1.4°  as  in  the  case  of  the  null  referen*  «,»  array.  The 
capture  effect  array  parameters  are  summarized  below. 


Lower  Antenna 

Middle  Antenna 

Height  * h 

Height  * 2h 

I - 1 

I„  » *0.5 

c 

c 

1 1 50  “ 0,34 

*150  " '°'08 

I90  - 0.46 

I90  - *0.32 

Upper  Antenna 


Height  * 3h 


I 

I 

I 


■ 

c 

isc 

90 


o 


-0.06 

0.06 
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Note  that  the  I and  Igg  current  amplitudes  given  above  represent 
the  suas  of  the  carrier  sideband  and  separate  sideband  signals. 
Again,  all  values  have  been  normalized  relative  to  the  carrier 
aaplitude  in  the  lower  antenna. 

It  should  be  noted  here  that  our  coaputer  prograa  automatic- 
ally offsets  the  eleaents  of  each  glide  slope  array  to  correct  for 
the  effects  of  proxiaity  phase  lag.  The  y coordinates  of  the 
eleaents  are  adjusted  so  that  all  dipoles  are  at  the  same  slant 
distance  froa  the  touchdown  point  on  the  runway  directly  opposite 
the  array.  For  each  array,  one  element  is  held  fixed  while  the 
other  eleaents  are  offset  relative  to  the  fixed  one.  For  the  null 
reference  and  sideband  reference  arrays,  the  lower  element  is  held 
fixed  while  fur  the  capture  array,  the  Biddle  element  is  held  fixed 
The  aaount  of  the  offset  for  each  dipole  can  be  calculated  approx- 
imately in  the  following  manner.  Let  y and  h denote,  respectively 
the  y-coordinate  and  height  of  the  fixed  antenna  element.  The 
y -displacement,  e,  of  any  other  element  in  the  array  relative  to 
the  fixed  one  is  given  approximately  by: 


e 


h2-H2 

~W7 


(66) 


where  H is  the  height  of  the  eleaent  which  is  to  be  offset. 

For  the  pi«;sent  study,  the  wavelength  X was  set  at  3 feet. 

All  arrays  were  positioned  300  feet  from  the  centerline  of  the 

runway.  The  heights  of  the  array  elements  were  set  at  14.33  feet 

and  28.66  feet  (approximately  5X  and  10X)  for  the  null  refeience 

* 

array,  7.17  feet  and  21.5  feet  (approximately  2.SA  and  7.5A)  for 
the  sideband  reference  array,  and  1.4.33  feet,  28.66  feet,  and 
42.99  feet  (approximately  SA,  iOA,  and  15 A ) for  the  capture  effect 
array.  Under  normal  circumstances,  a glide  angle  of  3°  would  be 
produced  by  all  three  arrays. 
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To  illustrate  the  performance  of  our  multipath  nodel,  the 
glide  slope  course  deragation  caused  by  a large,  flat  reflecting 
surfaee  (perhaps  the  side  of  a large  hangar  or  the  side  of  a hill) 
has  been  calculated.  The  structure  was  taken  to  be  300  feet  long 
and  100  feet  high  and  was  oriented  parallel  to  the  runway  center* 
line.  The  structure  was  positioned  200  feet  from  the  runway  center- 
line  in  the  y-direction,  on  the  opoosite  side  of  the  centerline  from 
the  glide  slope  array.  The  near  edge  of  the  structure  was  displaced 
1000  feet  in  the  x-direction  from  the  touchdown  point.  The  geome- 
try of  the  configuration  is  depicted  in  Figure  4. 

In  Figures  S,  6,  and  7 are  plotted  the  C.D.I.  signal  as  a 
function  of  the  horizontal  distance  from  touchdown  for  each  of  the 
three  arrays  for  the  arrangement  depicted  in  Figure  4.  In  each 
case,  the  receiver  followed  the  nominal  3°  zero  C.D.I.  hyperbola. 

The  results  for  all  three  arrays  are  very  similar.  Our  model  pre- 
dicts large  high  frequency,  alternating  fly  up  and  fly  down  signals 
near  the  end  of  the  flight  path  roughly  between  1700  feet  and  1800 
feet  from  touchdown.  This  region  in  fact  corresponds  almost 
exactly  to  the  region  of  specular  reflection  as  predicted  by  geomet- 
rical optics.  Note  that  the  excursions  are  large  enough  to  satur- 
ate the  receiver  over  much  of  the  affected  range  (full  s ale  de- 
flection * 150  microamps).  The  excursions  are  clearly  highly  localized 
and  damp  out  very  rapidly  as  the  aircraft  leaves  the  specular 
reflection  zone.  The  narrowness  of  the  pattern  (i.e.,  the  lack  of 
broadening  due  to  diffraction  and  the  rapid  damping  of  the  side 
lobes)  is  undoubtedly  attributable  to  the  large  size  of  the 
scatterer  relative  to  the  scattering  wavelength. 


GLIDE  SLOPE 


FIGURE  4.  OVERVIEW  OF  A MULTIPATH  CONFIGURATION 
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FIGURE 
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FIGURE  6.  SIDEBAND  REFERENCE  ARRAY,  FLYABILITY  RUN 
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FIGURE  7.  CAPTURE  EFFECT  ARRAY,  FLYABILITY  RUN 
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